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i-Q I We start from a field-theoretical model of zero range approximation to 

Q_i' derive three-dimensional relativistic two- and three-body equations on a null 



plane. We generalize those equations to finite range interactions. We pro- 
pose a three-body null-plane equation whose form is different from the one 



X 
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presented earlier in the framework of light-cone dynamics. We discuss the 

choices of the kernels in two- and three-body cases and apply our model to 
the description of meson and baryon Regge trajectories. Our approach over- 
comes some theoretical and phenomenological difficulties met in preceding 
relativized treatments of the three-body problem. 
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I. INTRODUCTION 

It is very well known that nonrelativistic potential models can be successfully applied 
to the description of the heavy qq spectra (see ref. [[jj] and the recent review ||). However 
when the states containing light quarks (u,d,s) are considered relativistic effects have to be 
taken into account ( see, e.g. |0-§[). Relativistic few body problem has received a great 
attention in hadronic and nuclear physics and a lot of papers have been devoted to the 
problem ( see e.g. reviews 0-0 and references therein). Different methods for deriving 
relativistic few-body three-dimensional equations have been discussed. Some authors use a 
diagrammatic approach, i. e., they select some leading diagrams and project them onto the 
three-dimensional momentum space. Others make use of effective Hamiltonians, employing 
various assumptions for their constructions. 

The most appealing are the relativistic approaches based on the null-plane (or light-cone) 
dynamics or, equivalently, on the analysis of Feynman diagrams in the infinite momentum 



frame (IMF) ( see recent papers ||T0| , |iTf and references therein). One of the main advantages 
of these approaches is the very well known fact that the wave function of the bound state 
has a simplest form in IMF, where pair creation from vacuum is suppressed. This was 
stressed a long time ago by Weinberg [Q. Indeed in the covariant representation Lorentz 



invariance of the S matrix at any order of perturbation theory is ensured by the contributions 
containing different numbers of particles and antiparticles in intermediate states. At the 
same time diagrams of noncovariant perturbation theory in IMF are dominated by the 
simplest intermediate states where all the particles have positive energies. 

To make calculations realistic we must restrict the number of degrees of freedom. This 
is usually done by imposing cutoff on the number of particles and considering only qq sector 



in the case of mesons and 3g sector in the case of baryons. In order to derive a two-body 
null plane equation in the three-dimensional form, one can use a Hamiltonian approach in 
which the null planes {x_ = -m (x° + x 3 ) = const} play the role of equal time surfaces 
J7, 10,11,13], or eliminate relative null plane time X- in the covariant approaches |7]||,|15|-17 . 



However in such approaches it is not trivial to define the angular momentum operators 
||7|j9|Jl3l| . This problem is circumvented in the papers which use the Relativistic Hamiltonian 
Dynamics (RHD) in the null plane form (or Null Plane Dynamics, NPD, or even Light Cone 
Dynamics, LCD), where the Poincare generators for the system of two or three interacting 
particles are directly constructed in terms of internal variables describing a system (see 

M). 

Generally there is no exact correspondence between the results based on RHD or NPD 
and quantum field theory (QFT). Such a correspondence can be found in the Zero Range Ap- 
proximation (ZRA), where the relative time can be explicitly eliminated and two-body three- 
dimensional null-plane equation can be derived from the Bethe-Salpeter equation p,|l7|. In 



this paper we discuss a generalisation of this equation to finite range interactions, which is 
equivalent to a representation of NPD, where the angular momentum operators have the 
same form as for free particles |13|], so that, under rather general assumptions, covariance 
of the theory is guaranteed. Although our results on the two-body null plane equation 
are not new, the discussion we present here gives new insight on the relation between NPD 
and QFT. Moreover we relate the kernel in the relativistic two-body null plane equation to 
the potential used by Godfrey and Isgur ||; this, together with linearity of meson Regge 
trajectories, suggests that for large separations between constituents the interaction should 
be of the oscillator type, whose parameter is proportional to the string tension squared. 



Such a kernel leads to predictions on Regge trajectories for mesons composed of light (u, 
d and s) quarks which are in good agreement with data, including small deviations from 
linearity for strange mesons. 

Another new point in our paper is the derivation of three-body null plane equations. 
Starting from the covariant three-body equation, we show that in ZRA two relative times 
can be eliminated and a three-dimensional three-body null plane equation for bound or 
scattering cases can be derived. We consider also a simple way of generalizing our treatment 
to finite range interactions, in line with Relativistic Hamiltonian Dynamics. 

The form of the three-body null plane equations is different from the one proposed 



earlier in the framework of RHD |fL8 ], as well as from the naive relativistic generalization 



of the Schrodinger equation used in refs. [|3],[|,[19| for the description of baryon spectra. In 



particular, as we show in the present paper, our approach presents theoretical and also 
phenomenological advantages over these treatments. 

We apply three-body null plane equations to Regge trajectories for a system of three 
relativistic quarks, analyzing the kernel corresponding to a three-body force with string 
junction, involved in the mass operator squared. We take this operator in a form which 
commutes with angular momentum in a representation which coincides with the one for free 
particles, so that also the three body equation yields a covariant [|7| description of dynamics. 
Approximate linearity of Regge trajectories fixes the dependence of the kernel on two differ- 
ent relative coordinates, which for large relative separations has to be still of the oscillator 
type. Equating the slopes of two Regge trajectories corresponding to the orbital excitations 
of two relative degrees of freedom fixes the sector that describes the relativistic recoil of the 
moving two-quark subsystem. The implementation of the condition that the slope of the 



diquark-quark orbital excitations should be same as for meson Regge trajectories fixes the 
oscillator parameter. We apply the model to the baryon Regge trajectories composed from 
the light quarks and show that the predictions of our model are in good agreement with 
data, reproducing small deviations from linearity. 

It is worth stressing that, unlike the approaches based on a relativized Schrodinger equa- 
tion, our three-body Hamiltonian satisfies cluster separability for two-body forces; further- 
more, implementing our model with three-body interactions allows to pass from mesons 
to baryons, without re-adjusting fundamental parameters or assuming a different kind of 
interaction between quarks. 

The paper is organized as follows. In Sect. 2 we carry on the program of studying the 
relativistic two-body equation on a null plane. Starting from the covariant Bethe-Salpeter 
equation and from a field theory model in ZRA, we derive a relativistic two-body equation 
for a null-plane wave-function. Then we generalize this equation for finite range interactions 
and introduce the relativistic null-plane Lippmann- Schwinger equation. 

In Sect. 3 we focus on the three-body equation. As in the two-body case, we start from 
ZRA; then, using the form of the equation introduced for finite range interaction in the 
two-body case, we formulate a relativistic equation for a three-body system appropriate for 
this kind of interactions. Moreover we show that this equation satisfies cluster separability 
condition. 

In Sects. 4 and 5 we apply two- and three-body equations to the study of meson and 
baryon Regge trajectories. We discuss the choice of the kernels in two- and three-body 
cases and we apply the model to the description of the experimental data of meson Regge 
trajectories composed of qq, qs (or sq) and ss orbital excitations and Regge trajectories of 



baryons composed of u- and/or d-quarks. We consider also the A Regge trajectory. 
Sect. 6 is devoted to conclusions. 

II. TWO - BODY NULL PLANE EQUATION 

Let us consider the case of two non-identical, spinless interacting particles, with masses 
m\ and m 2 respectively. The Bethe-Salpeter (BS) equation for the reduced amplitude \ 



reads 20,21 



(p? - ml)(pl - ml) X {pi,P2) + J 0^^i^)^ P ~ Pi ~ P2)Ui2 (P,P'; P)X(P' 1 ,P' 2 ) = 0, 

(2.1) 



where 



P = VPi - (1 - V)P2, P' = VP'i - (1 - V)P2 (° < V < 1) ( 2 - 2 ) 

and U12 is the kernel in momentum space. 

The BS equation meets some intrinsic difficulties. First of all, the normalization condition 
for x 

(2.3) 

depends on the interaction through the kernel Ui2(p,p';P), which is a source of complica- 
tions. 

Moreover x depends on the relative energy, therefore Xp{ x ) depends on the relative time 
t. This property of the reduced amplitude x leads to different complications when we try 
to apply the BS equation to physical problems. Indeed, in the ladder approximation the 
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BS approach does not exhibit the correct static limit, i. e., when, e. g., mi becomes large, 
while keeping 777,2 fixed: some solutions exist with negative norm or which do not satisfy the 
correct non-relativistic limit. This problem can be avoided by eliminating the relative time 
dependence PJT?]. 



Let us define a null plane projection of the BS amplitude as follows: 

"0(Pi,P2) = / dpi^dp 2 -Pi + p 2 +S(P- ~P\- -P2-)x(Pi,P2) , (2.4) 

where p, are null plane three-momenta, i. e., 

Pi = (PiT,Pi+) (2.5) 



and 



Pi+ = J=—, PiT = {PtxiPiy) • (2.6) 



We assume the kernel of eq. fl2.1|) to be independent of relative null-plane energy p_, i. 



U 12 = U 12 (p,p';P), (2.7) 



where 



P = Wi - (1-^)P2, p' = ^p'x - (1 - 77)p' 2 . (2.8) 

This amounts to assuming an instantaneous interaction in the null plane time x+. Then the 
reduced BS amplitude may be rewritten as 

, , l /rfr 12 [/ 12 (p,p';P)^(p'„p' 2 ) 

* { ?^ = - W-^Xrf-m?) ' (2 ' 9) 

where 



dT 12 = ^^^ 3 (P - P \ - p' 2 ). (2.10) 

2(27r) 3 p 1+ p 2+ 



Integrating both sides of (|2.9| ) over Pi+p 2 +dpi_dp 2 _5(P- — p\- — P2-), we derive the 
following equation for a null plane wave-function: 

V>(Pi> P2) = iP+G (xi, Pit; x 2 , p 2 r) y dr 12 £/ 12 (p, p'; P)V>(p'i, P 2 ), ( 2 - n ) 

where 

G (*i, Pir; x 2 , P2T ) = -[PktM + Pk±^l _ M 2 - P T - .e]- 1 (2.12) 

Xi x 2 

and 

M 2 = P 2 , x« = ^, Pt = Pit + P 2 t- (2.13) 

Moreover let us assume the kernel U\ 2 to be of the type 

U 12 = U 12 {x,p T ;x',p' T ) , (2.14) 

where 

a; = tjxx - (1 - r])x 2 , Pt = VPit ~ (I ~ v)P2T ■ (2.15) 



Then we may define a null-plane wave-function \l/ that, according to fl2.11|) , satisfies the 
following equation: 



" ,2 + Pir + ^l _ M ,..p, 1xRt . |-pr:-< . 2p ., y , ) ., 



rPlT + m l , P 2 T + m 2 i\j2 T >? 



Xi x 2 



i J dL' 12 U 12 (x, p T ; x', p' T )^(x[, p' ir ; a£, p' 2T ) = °; ( 2 - 16 ) 



where 



1 2 r/r' 

dL '™ = 979^tn -T^t]^(1 - *i - x' 2 )5 2 (P T - p' 1T - p' 2T ). (2.17) 



At this point we discuss approximation (|2.7Q . The simplest case is ZRA, i. e., when U12 
is a constant. In this case the null plane wave-function reads 

N 
* = — g 2 2 > ( 2 - 18 ) 

where we have introduced the variables 

£ = xi, qr = (1 - OPir - ^P2T (2.19) 

and the constant N is determined by the normalization condition ( |2.3| ), i. e., 

w/^ |WK '* )|, = 1 " (2 - 20) 

which turns out to coincide with the usual normalization of the infinite momentum frame 
wave-function. 

In this limit a null plane description is equivalent to field theory; however it has limited 
applications. It is therefore important to construct kernels of finite size, using different 
models. For example, it is possible to use ladder approximation to the BS equation and to 



project it on the null plane. Such an equation was suggested by Weinberg |TJ]. Another 
example was considered by t'Hooft |L5]], who derived a two-body null plane equation for the 
qq system in 1 + 1 QCD, taking the limit of N c — > 00, in such a way that a s N c have a finite, 
constant value. 

A great deal of papers were devoted to the analysis of different aspects of Hamiltonian 
null plane approach flnj. Here we try to impose the Tamm-Dancoff cutoff on the number 
of particles, while conserving a very simple representation of angular momentum operators. 
Generally speaking, in null plane formalism angular momentum operators should depend 
on interaction @H; in particular this is the case of the Weinberg equation J12J , where the 



representation of angular momentum operators cannot be defined straightforwardly, since 
these operators would depend on interaction in a very complicated way. Furthermore, as 
stressed by Leutwyler and Stern [7| , without defining angular momentum, the wave function 
cannot be normalized in a usual way. Therefore the normalization condition for the wave- 
function satisfying the Weinberg equation is not clear, except in ZRA, where the S-wave 
alone is involved. 



However, as shown by Terent'ev ||13|| , in the two-body sector with a finite range of 
interaction it is also possible to choose for the angular momentum operators the same rep- 
resentation as for free particles. For example, in the relativistic oscillator model |2I1 the two 
different representations are related by a "gauge" transformation ||. Usually one assumes 
this equivalence to hold true also for other interactions. In the Terent'ev representation the 
angular momentum operator takes the form 

f=-*qx ^, (2.21) 

aq 

where 

q=(qT,g,), ^ = ^-^-^p ( 2 -22) 

M 12 = M 12 (q) = ei + e 2 , (2.23) 



cp + mf, z = l,2 (2.24) 

and the variables £ and qx are defined by eqs. ( |2.19| ). In this representation angular mo- 
mentum conservation is ensured by rotation invariance of the kernel in eq. (|2.16|) . 
Then it is possible to rewrite eq. ( [2.16|) in the form 



Ml 2 m = M 2 ^ , (2.25) 
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where 



M? 2 = MUq) + W 12 , 



(2.26) 



and 



W-12 = iUio. 



(2.27) 



If W12 commutes with the angular momentum operator, the angular condition is satisfied 
and a covariant |7| description of the composite system is guaranteed. 

Extension to the case or two interacting spinning particles can be performed in a straight- 



forward way by using the Terent'ev representation [13]. In this case the meson state vector 
can be written in the form 



P+i Pt, J, A >np- 



E E C\l Xa CYr f^%(q)Y ln (q)%,,(q L )t' V (-q)\pifi',P-piU > > NP , 



where NP means "null plane", the C's are Clebsch-Gordan coefficients, Y[ n the Legendre 
functions, q is defined by eq. ( |2.22| ), pi = pi(q) is the four-momentum of the quark, 



q=\q 



« q 

q = - 

q 



(2.28) 



and the V^'s are the Melosh [^] rotation matrices for quark and antiquark, i. e. |L3 



V(q) 



V(-q) 



m 1 + e 1 + q z + ie js (jjq t 
[2(e 1 +m 1 )(e 1 + q z )]* 



m 2 + e 2 -q z - ie js ajq s 



(2.29) 



(2.30) 



[2(e 2 + m 2 )(e 2 - q z )\* 

The state vector | P + ,Pt, J, A >np at four-momentum P is obtained from the state 
vector at rest | P 0+ , Ot, J, A >np by the element of the Poincare group l(P <— P ) equivalent 
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to the product of two Lorentz transformations, i. e., X(P <— P^) X(Poo <— Pq), where 
P = (M, 0,0,0). 

Now let us consider nucleons interacting through exchange of different mesons. The 
kernel that describes such an interaction can be identified with the Born term expressed 
in terms of null - plane variables. For example, in the cases of scalar and vector meson 



exchange, the kernels have, respectively, the following forms ||14|| : 



W 12 ., V = V + $)WV(Z', q)V(q), 



(2.31) 



where 



V(q) = V( Pl ) ® V(p 2 ) , V(q') = V(p\) ® V(p' 2 ) 



(2.32) 



^S(q',q) 



2 i ,.2 



(q'~ q)' + /i 



u(p[)u( Pl ) u(p' 2 )u(p 2 ), pP = (ei' } , (-) 1+i q (,) ), (2.33) 



r{r) 



(Wl 2 i)a>f3>,ap(<l', q) 



(q'- q) +/i 



o 2 ^(q'. q)K^(-q' ; -q) sv> 



(2.34) 



Vlj(q', q) = (q' + q)<W + i (<J - q) X a . (i = 1, 2) 



a/3 



(2.35) 



Oq',q) 



v e + m) + q ■ q' 



<W + ia a p • q x q' 



e + m 



(2.36) 



fj, is the mass of the exchanged meson and q' is defined analogously to q (eq. [2.22|) , 
with Myi substituted by M[ 2 , which generally assumes a different value. In the case of 
vector meson exchange we have assumed the two interacting particles to have equal masses, 
rrii — m 2 — m i therefore we have also C\ = t 2 = e. In the Terent'ev representation of angular 



momentum, such kernels satisfy the angular condition |18 |. 
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If we consider kernels of this kind, there is no exact correspondence between null plane 
dynamics in two - body sector and field theory; however the former description may be a 
good approximation to the latter one if only two - body intermediate states are important. 
A further improvement would consist in constructing multichannel null - plane dynamics; 
indeed, if the number of channels becomes infinite, null - plane dynamics and field theory 
will turn out to be equivalent. 

Equation ( |2.16|) may be also written symbolically as 



* = G W 12 ^, (2.37) 



where Gq is defined by eq. (|2.12 ). This leads immediately to the Lippmann-Schwinger 



equation for the two body relativistic T - matrix, i. e., 

T = W 12 + W 12 G T, (2.38) 

which, in momentum representation, reads 

T(xi, Pit, x 2 , p 2 r, x[, p' 1T , x 2 , p' 2T ) = W 12 (x, p T ; a/, p' T ) + 

f, r „ W 12 (x, p T , x", p" T )T(x", p" 1T , a?2, p" 2 ; x\, p' 1T , x' 2 , p' 2T ) 

' / 12 „ n n" 2 4-m 2 n" 2 -4-m 2 ' {^.Oj) 



As proved in ref. 13H , in this scheme angular momentum operators are the same as for 



free particles, the only difference consisting in the appearence of Melosh rotation matrices 
in partial wave expansion of the wave-function and of the scattering amplitude. 

III. THREE-BODY NULL PLANE EQUATION 

Analogously to the two-body case, the reduced three-body amplitude reads 

13 



(Pi - m l)(P2 ~ m 2)(P3 ~ m -i) ( f ) (PliP2,P3) + 



+ ' ^ t^Y ^Y (2?r) ^ 4(P ~ P[ ~ A ~ ^ U ™^ n " P )^l' *^ = ' ( 31 ) 



where C/103 is the kernel and 



n = (pi , p 2 , ps) , n = (pi , p 2 , p 3 ) 



(3.2) 



The null plane wave-function is defined by 



^'(Pi>P2,ps) = / dp 1 „dp 2 -dp 3 ^pi + p 2+ P3+5(P- -pi_ -p 2 - -p 3 -)0(pi,p 2 ,P3) • (3.3) 



Now let us suppose the kernel to be independent of p r _, i. e., 



£/ 12 3 = £/i23(n,n' ; p), 



(3.4) 



where 



n = (Pi,p 2 ,p 3 ), n' = (p'l.p'a.p's)- 



(3.5) 



In this case we derive from (|3.1| ) the following equation for null plane three-body wave- 
function: 



V'(II) = -G'^pir) J dL' l23 U 123 (II,II>;P)i>>(Il> 



(3.6) 



where 



3 A^l 



dx' 
dL 'vn = A(0 n 6 II — r d P'rT<*(l - x'x - x 2 - x' 3 )S 2 (P T - Pit - P 2 t - Pst) (3-7) 



and 



G (Xi, PiT) 



PlT + m l , Plr + m 2 , V>l T + m l 



■I'l 



+ 



;r 2 



+ 



M 2 



Pl-ie 



■1-2 



(3.8) 
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If we assume the kernel U123 to depend on p r+ only through fractional momenta x r , the 
three-body equation reads 



Pit + m i 1 P 2 r + m 2 , P 3 r + m 3 _ m 2 — P 

X\ ' X2 X3 



^'(^,PiT) + 



+ / dL' 123 U 123 (x r , p rT ; a£, p'^^'K, P'oO = ° 



(3.9) 



We write, as in the two-body case, the Lippmann-Schwinger equation for three body T- 



matrix, 1. e., 



T(xi, p lT ] x'i, v'iT) = U 12 3(x r , p rT ; x' s , p' sT ) + 



+ dL" 



123 



M Z + P Z T - 



U\Zz{x r , PrT, x' s} p sT )T(x" , p" iT ] X-, p iT ) 

2 p"\ T + m\ p"\ T + m\ p"l T + m 2 3 



1 1 



x 2 



1 3 



+ ie 



(3.10) 



Now we assume only two-body interactions to occur. Then the kernel of eq. (|3.9|) may 



be written as 

Ui23{Xr,Prr,x' s ,p' sT ) 



i+j+k 



XkUijixy, p.. T ; x' tj , p' )2(2ir) 3 5(x k - x k )5 z {p kT - p' kT ) 



'ijT^ijir ijT 



(3.11) 



where p^r and Xij are defined analogously to eq. Q2.15Q . Then equation (|3~9| ) can be written 



as 



Mj^tf' = MV, 



(3.12) 



where we have set 






+ .£. T 



Wi. 



ift¥=k x Xfc 



Wij*' = i / dL , l' j U lj (x ij1 p ijT ; x' i:j , p' ijT )^'(x[, p' iT ; x'j, p' jT ; x k , p kT ) 



(3.13) 
(3.14) 
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and dLfj is denned by (|2.17| ), substituting the overall transverse momentum P T by 

PijT = PiT + PjT- (3.15) 



In the case of ZRA W^ are constants and condition ( |3.4| ) is satisfied automatically. Such 
an equation is equally satisfied if the two-body interactions are of the type discussed in the 
previous section. Here we show that eq. (|3.9| ), with a kernel like (|3.11|) , satisfies cluster 



separability condition. Indeed, if, e. g., the third particle is taken far away from the other 
two, we have W 2 j, = W31 = and eq. ( ^.9| ) reduces to 

P 2 12T + M? 2 + p 2 3T + ml _ = M 

1 - x 3 x 3 

where M\ 2 * s given by eq. (|2.26| ). 

Dividing both sides of eq. ( p,16|) by 2P + , we get 

P_*' = (P 12 _+p 3 -)*', (3.17) 

which proves cluster separability for P_. 

In this connection let us note that the semi-relativistic three-body Hamiltonian 

H = T(pi, p 2 , p 3 ) + V 12 3 (pi, p 2 , P3) , (3.18) 



where T is the relativistic kinetic term, 



T(Pi,P2,P3) = £ [p-+rr^Y (3.19) 

t=i 

and the potential is assumed to be the sum of two-body interactions, i. e., 



Vl23 = ZX ( 3 - 20 ) 

i<j 
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does not satisfy the cluster separability condition. Indeed, if one particle is taken far away, 
the Hamiltonian in the centre-of-mass system should have the form 



Tjc.m. 
-"12+3 



Q?2 + M 1 2 2 + V /Q? 2 + m|, 



(3.21) 



where Q12 is the total momentum of the pair (1 — 2). 



On the other hand, in ref. |18| the three-body mass operator has been assumed to be of 



the form 



TTC.m. 

-"123 



M + J2(Eij - Ey 



(3.22) 



i<j 



where Mq is the overall free mass operator, 



E lj = x Ql + M^ j E lJ = J(^ 1 + M i 



iij 1 x "ijj 



(3.23) 



and M?j, M 2 - are defined analogously to eqs. ( 2.26 ) and ( |2.23 ). Evidently the mass 
operator ( |3.22|) satisfies the cluster separability condition ( p. 17] ). However, if we define the 
mass operator squared for the three-body system as 



M 10 q = (H^-) , 



'123 



(3.24) 



this is different from eq. ( [3.13J ), derived from the BS equation. 

It is worth noting that the mass operator ( |3.22j ), which was considered in |n|, commutes 
with the angular momentum operators taken in the same representation as for free particles. 



This is not true for the mass operator defined by eq. ( |3.13|) . However we do not worry about 
this; indeed, when considering the baryon Regge trajectories in Sect. 5, we shall propose a 
different kernel, which correponds to a three-force with string junction and commutes with 
the free angular momentum operators. 

We note also that eq. ( |3.13| ) has been considered in ref. [23], however the form of this 
equation has been postulated, no argument has been given for deriving it. 
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IV. MESON REGGE TRAJECTORIES 

Now we apply two- and three-body equations to the description of meson and baryon 
Regge trajectories. In particular in this section we consider mesons. We relate the qq 
interaction term Wi 2 (eq. Q2.26| )) in the two-body equation to the potential used by Godfrey 



and Isgur [RL i. e., 



V l2 = V™ nf + V$ p + V°° (4.1) 



where 



-rrCOn/ 

V X2 



4 a„(r) 

c + vr 



3 r 
c = -0.253 GeV, v = 0.18 GeV\ (4.2) 

r is the distance between q and q in the two-body cms and V r 12 w and V{£ are, respectively, 
the hyperfine and spin-orbit interactions 0. Squaring the mass operator 

M12 = M 12 + V12, (4.3) 

we recover eq. Q2.26J ), with 

W 12 = M 12 V 12 + V 12 M 12 + Vl 2 . (4.4) 

Obviously in this case eq. ( |2.26|) provides the same spectrum as Godfrey and Isgur 0. 

Here we concentrate our attention on Regge trajectories for systems of light quarks. 
For large quark- ant iquark separations the hyperfine and spin-orbit interactions (except for 
the Thomas precession term) can be neglected and the main contribution comes from the 
confining term. On the other hand, we know from experiment that Regge trajectories are to 



a high extent linear, with the same slope. For large orbital excitations a linear /-dependence 
(where I is the orbital angular momentum) of squared masses, i. e., 

M 2 = 8z/Z (4.5) 

can be found by using a linear confining potential in the mass operator. Now we pose 
the question how to reproduce phenomenologically such a linearity in the framework of the 
two-body equation (|2.25|) . If we assume Wu to be of the form 



W 12 (r) = wV, (4.6) 

then for massless constituents and for orbital angular momenta I » 1 we may consider the 
semiclassical limit for M 2 , i. e., 
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M 2 = 4— + cA c . (4.7) 

The mass squared of the system has a minimum when r satisfies the equation 

dM2 0, (4.8) 



dr 
whose solution is 



In this case we have the following Regge trajectory: 



r = (-^)— 2 - (4.9) 



or Qui 

which turns out to be linear when ( = 2. In this case 

M 2 = Aul. (4.11) 
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Comparing eq. ( f4.11|) with eq. (|4.5|) , we find uo = 2u. 



This means that for large q — q separations the kernel W\2 is of the oscillator form. 
Taking into account that all the Regge trajectories are approximately linear, with the same 



slope, even for low I, (see, for example, [26,27]), we can set the kernel A4.4Q in the form 



W 12 = W + 4z/ V + AW 12 (r) , (4. 12) 

where AWulr) - which generates non-linear corrections to Regge trajectories, e. g., spin 
dependent terms - is expected to be small. 

We apply eq. ( |2.25| ) with the kernel ( |4.12| ) to the description of qq Regge trajectories, 



where q = u, d or s. As it follows from the calculations of qq meson spectrum done by 
Godfrey and Isgur ||, spin-dependent corrections are quite important for Is and 2s-states. 
For lp-states they should also be taken into account: the splitting between 1 3 P 2 and l 3 Po 
levels is about 200 MeV. However for all other states those corrections are small and can 
be neglected (see also the discussion at the end of this section). If we neglect AWi 2 (r), the 
spectra of mesons composed of u and d quarks are given by 

M 2 Z = Am 2 + W + Au[2n + I -1/2] (n > 1) . (4.13) 

We take the same masses of u and d (constituent) quarks as in ref. 0, i.e. , 

m = m u = m d = .22 GeV. (4.14) 

The two parameters Wq and to have been fixed fitting the p-trajectory (see Fig. 1(a)), i.e., 

W = -1.255 GeV 2 , u = 0.2836 GeV 2 . (4.15) 

There are eight meson Regge trajectories made of u and d quarks. They are degenerate in 
isospin. For each isospin / there are three trajectories with angular momenta j = I + 1, 
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j = I and j = I — 1 and total quark-antiquark spin su = 1 and one trajectory with j = I 
and S12 = 0. We use the notation 

M/ =J ._ 1)Si2=1 , M[ =jtna=1 , M[ =j+lt8la=1 , Ml hSi2=Q . (4.16) 



All such trajectories are plotted in Figs. 1. Solid lines are determined using eq. ( |4 . 1 3[ ) . We 



observe that all available experimental data are in good agreement with the spectrum given 
by eq. (|4.13|) for / > 1. Approximate linearity with equal slopes of all nonstrange meson 
Regge trajectories for I > 1 confirms a posteriori that nonlinear, spin dependent terms of 
the potential fall down rapidly at increasing I. At I = these contributions have to be 
taken into account, as appears in the plots. Such deviations are smallest in the p-trajectory, 
suggesting that in this case nonlinear effects are less important. 

When the masses of the quark and antiquark are different, we represent the mass operator 
squared in the form 

M\ 2 = 4(aq 2 + /i 2 ) + W 12 + G{a, fi 2 , q 2 ) , (4.17) 



where 



rri] + mo 
/i= 1 2 (4.18) 



and 

G{a, fi 2 , q 2 ) = (y/q 2 + m 2 + y/q 2 + m 2 ) 2 - 4(aq 2 + /i 2 ) . (4.19) 

Note that G is equal to when 

• m\ = rri2 and a = 1 ; 

• q — > 0, for any value of mi and m 2 ; 
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q > > m\ , m 2 , and a = 1 . 



If mi 7^ m 2 , the eigenvalues and eigenfunctions of the operator ( |4.17| ) can be suitably 
approximated by solving the equation 



4 (aq 2 + /i 2 ) + Wvi 



m = M A m 



(4.20) 



where W\z is given by eq. ( |4.12| ), with Ajy 12 = 0. We choose a in such a way that G 
corrections to eigenvalues of ( |4.20| ) vanish to first order. This can be realized to a good 
approximation by imposing a to be a solution of the equation 



G(a, [i , q 2 (a)) = 



(4.21) 



where q 2 (a) is defined as the average value of q 2 over the wave function \1/, which obviously 
depends itself from a. As can be checked, eq. ( |4.21|) is an algebraic equation of the third 
degree in yfa\ among the roots of this equation we pick up the one which tends to 1 as 
m\ — m2 tends to 0. Let us examine the solution in detail. The eigenfunctions of eq. (f4.20|) 



are 



tfnZm(r) 






G 



ni — C n ie 

fjl+3/2 



l/2/3V r ,y 



(/3r)'F(l-n,J + -,/3V 



ni 



r(z + 3/2)\ r 



2T(n + I + 1/2) 



n 



d 2 



U! 



2v^' 



q 2 (a) depends on n, I: 



q 2 (a) =< nlmlq^lnlm >= (n + 1/2 — 1/4) —= 



(4.22) 
(4.23) 

(4.24) 
(4.25) 



(4.26) 
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having denoted by \nlrn > the eigenstates of eq. ( |4.20| ). Therefore the solution of eq. ( f4.21|) 



which we have picked up also depends on n, I through the combination g n \ = In + / — 1/2. 
Then the eigenvalues of equation (f4.17|) can be approximated by the following expression: 



M l nl = V + W + 4 (2n + I - 1/2) cu^i . (4.27) 

The plot of a n i as a function of g n \ is shown in Fig. 2 for m\ = .22 GeV and irt2 ranging 
from .22 to 1.22 GeV. When m 2 < .7 GeV, the difference between a n i and 1 does not exceed 
10%. In this case the correction to the slope of the Regge trajectory does not exceed 5% 
and the deviations from linearity are negligible. 

We apply (|4.27| ) to the description of the strange meson trajectories. For the mass of the 



strange quark we used the same value as in ref. ||, i.e. m s = .419 GeV. Then we can predict 
all the strange meson trajectories {M l ' ) and mesons with hidden strangeness (M? ) - see 
Figs. 3 and 4. 

As in the case of non-strange mesons, the agreement between our theoretical predictions 
and the existing experimental data |2S| is good, except for the states with / = 0. 



The spectra given by eq. ( 4.27 ) contain also radial excitations. The mass of the qq state 
with n = 2, I = is predicted around 1.7 GeV, while the in pseudoscalar channel the first 
radial excitation is around 1.3 GeV (7r(1300) and 77(1295)) and in the vector channel it is 
around 1.4 GeV (u;(1390) and p(1450)). This discrepancy can be due to neglecting, in our 
model, spin-spin interaction, which mixes states with different n. 

V. BARYON REGGE TRAJECTORIES 

Let us apply now the three-body null plane equation ( |3.13| ) to the description of the 
baryon Regge trajectories. We consider a system of three light quarks in a colour singlet 
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state. For orbital excitations with large angular momenta we may distinguish among three 
different configurations: 

f . String-like configuration, with two quarks forming a diquark ( r\i ~ r 2 3 ^> r%i, where 
Tij is the relative distance between the i— th and the j— th quark in the two-particle 
cms) ["D" -configuration]; 

2. Symmetric triangle-like configuration ( ri 2 ~ r 2 3 ~ r^\ ) ["T" -configuration]; 

3. Star-like configuration with string junction [ "Y" -configuration] . 

In the first case, the problem reduces to the two-body case; in particular u is the same as 
in quark- ant iquark system, since in a baryon the diquark is in a colour anti-triplet state. For 
rather large overall angular momenta, when the rest mass of the diquark may be neglected, 
the Regge trajectory is the same as in two-body case, eq. ( |4. 1 1| ) . 



In the second case we recover eq. ( [3.131 ), where the two-body kernels Wij are of the 
oscillator form. Taking into account the colour degree of freedom yields 

Wij = \o\%. (5.1) 

If, moreover, we take into account of the symmetry of the configuration and proceed 
similarly to the two body case, we find 

M 2 = 6cuL 123 , (5.2) 

where L 123 is the overall angular momentum of the three-quark system. Therefore for a 
fixed L 123 the diquark-quark configuration has a lower energy than the symmetric triangle 
configuration, according to the following ratio: 
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M^:M 2 D = 1.5 : 1. (5.3) 

In the case of a relativized three-body Schrodinger equation, i. e., 



E Wi + m i + V ™)*' = M*', (5.4) 



»=i 



where V123 is linear either with respect to relative distances between quarks (triangle con- 
figuration) or with respect to distances of quarks from the string centre (string junction 



configuration), the situation is different; in this case it results |29 



Ml:M 2 D = 1.3 : 1.0. (5.5) 

The difference in the ratio M^ : M|, is due to the relativistic recoil effects taken into account 
in the three-body null plane equation. 

In connection with triangle configuration we recall, as shown in sect. 3, that the three- 
body semi-relativistic Hamiltonian with a two-body potential does not satisfy cluster sep- 
arability condition. Alternatively, using three-body mass operator (|3.22|) - which satisfies 



cluster separability condition - would lead to a non-linear Regge trajectory. In fact, if we 
assume the potential Vy, which describes quark-quark interaction in the two-quark cms, 
to be proportional to the distance between the two interacting particles - as imposed by 
linearity of meson Regge trajectories -, the same potential in the overall cms is no longer 
linear, since, according to eqs. ( |3.22| ) and ( |3.23| ), it transforms to 



VL = E ii - E v ( 5 - 6 ) 



where 



Eij = \/Q|- + {Ma + V %j f , Eq = \/Ql + M£ (5.7) 
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and Mij is the two-body free mass operator, defined analogously to (|2.23|) . 

In the case of a star-like configuration with string junction the interaction is not a two 
body one. This is a genuine three-body force. Indeed arguments in favour of a three-body 



force for describing the baryon spectrum have been illustrated in ref. [30]. Also in the 
framework of a relativized Schrodinger equation a three-body force with string junction ||] 
allows a unified treatment of mesons and baryons, whereas exclusive use of two-body forces 
H demands a quark-diquark structure. Then we assume a kernel of the form 

U 123 = U + u ,2 (r 2 01 + r 2 02 + r 2 03 ), r 0i = \r 0i \, (z = 1,2,3), (5.8) 

where Uq is a constant and ?oj is the vector from the string center to the i— th quark in the 
overall cms. It is not difficult to recognize that such a kernel is of the type ( |3.4| ), since r j 
is conjugate to momentum 

Poi = Pi - Po, (5.9) 

where 

Pi = (Pit,Pi z ) (/ = 0,1,2,3) (5.10) 

and the pi z are related to the fractional momenta £/ by means of the following equations: 



6= " m ' + M +PlZ , = 0,1,2,3). (5.11) 

The kernel ( |5.8| ) contains also the dependence on cm. coordinate. To eliminate it, we 
proceed as in the non relativistic case of three particles with equal masses, defining the 
relative coordinates as follows: 

R = g (roi + ro2 + ro 3 ) 
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p = r 01 - r 02 



(5.12) 



A = - (roi + ?"02 - 2ro 3 ) 



In order to take into account relativistic effects, we write the interaction term in the form 



-? / 2 r 2 . 1 -e 



#123 = U + J ( -?7aA + -7] p p 



(5.13) 



In the non relativistic limit the factors rjx and r\ p tend to 1. Let us stress a difference between 
our model and the non-relativistic harmonic oscillator model. In our case the interaction 



term ( {5.13Q appears in the mass operator squared, i. e., 



^123 



Q 2 + M 1 2 2 (q2) + v /Q 2 + m| 



+ U 



123, 



(5.14) 



(we have used shortened notation Q for Q12), while in the non relativistic case it is a part of 
the mass operator. Therefore the mass spectrum and the Regge trajectories in the relativistic 
and non relativistic models with harmonic oscillator interaction will be essentially different 
for systems with light quarks. The momenta Q and q in equation (|5.14j) are, respectively, 
the momenta of the subsystems (12) + 3 and 1 + 2 in their cm. frames. 

In the nonrelativistic case A and p are the conjugate coordinates, respectively, to the 
momenta Q and q. In the relativistic case this can be justified only for A. Of course the value 
of I A I is not Lorentz invariant and in the overall cm. system it can be regarded as the distance 
between the subsystem (1 + 2) and particle 3. In line with these considerations we can take 
the coefficient i]\ in eq. (|5.13|) equal to 1. The parameter r] p , which describes the relativistic 
recoil effect for subsystem 1 + 2, will be fixed later. Therefore in the relativistic case the 
choice of the kernel ( 5.14 ), where both A and p are considered as conjugate coordinates to the 
momenta Q and q should be treated as an " ansatz" , which however does not correspond to 
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a sum of two-body interactions. An important consequence of this choice of the interaction 
is that it commutes with the angular momentum operators taken in the same representation 



as for free particles: 



'123 



7 7 * d - d 

<9Q <9q 



(5.15) 



Since the mass and angular momentum operators have been fixed according to eqs. ( |5.13| ) 
and ( |5.15|) respectively, all other Poincare generators can be constructed according to the 
scheme prposed by Berestetsky and Terent'ev |H] (see also [0). The angular condition is 
satisfied, therefore the eigenstates of our Hamiltonian are guaranteed to be the null plane 
restrictions of covariant wavefunctions [|7| . 

Let us consider the spectrum predicted by the mass operator squared ( |5.14j) . Using the 
results of the previous section, we write the kinetic term in eq. ( p.!4| ) as follows: 



Q 2 + M 1 2 2 (q2) + v /Q 2 + m 



4(«Q 2 + /i 2 ) + G(«,/i 2 ,Q 2 ) 



/i = -[M 12 (q)+m 3 ] 



(5.16) 

(5.17) 



and we assume ni\ = rri2 = m. As in the two-body case, we approximate the eigenvalues 
and eigenfunctions of the operator M 2 23 ( eq. (|5.14f) ) with those of the eigenvalue equation 



v' ,,2 , , ,' 2 ( \2 , 1 ^e 



4(aQ 2 + /i 2 )+^ z (-A 2 + -^) ■■(/<, 



m = M 2 m 



(5.18) 



where a solves the equation 



G(a,^ 2 (q 2 ),Q 2 ) = 0. 



(5.19) 



Here q 2 and Q 2 are the average values of q 2 and Q 2 for the eigenfunctions of eq. ( |5.18 ) 
which can be written in the form 



2N 



^NLM,nlm (Q, qj = V'VLAf (Qj </>raim (q) , (N,n> 1) (5.20) 

where the functions ^nlm (Q) and n ; m (q) are defined by the expressions ( [L22H4.25| ), with 



1 J 



J 



^=7575- # = 2?5^ (5 ' 21) 



Assuming M 2 2 (q) >> m 2 , and moreover that [To = 2W , where Wo is the same as for gg 
system (which amounts to taking Wq per degree of freedom), we can write the eigenvalues 
of eq. (|5.18| ) as follows: 



M 2 NLnl = M 2 nl + W + 4u/y ~a NL (2N + L- 1/2) , (5.22) 

where 

M 2 Z = W + Am 2 + 2V2cu'^f p (2n + / - 1/2) . (5.23) 

Let us consider the limiting configuration of L — > oo and I fixed, which corresponds to a 
string-like diquark-quark configuration. In this case, since Q 2 >> M 2 2 ,ml and «tvl — > 1, 



the slope of the Regge trajectory for this configuration results to be */3/2(4o/) 1 ; this slope 
equals the slope of the meson trajectory, therefore 

u' = ^3/2uj . (5.24) 

On the other hand, in the limit I —>■ oo, with L fixed, the slope of the Regge trajectory results 
to be (2\/3u}^f}^)~ 1 . As the two slopes are equal, r\ p is uniquely fixed to 4/3. Therefore the 
mass spectrum ( 5.22| ) has the form 



M 2 N Lni = {-2W + 4m 2 + 4u[yfa{2N + L)+2n + l - 1/2(1 + yfie)]} GeV 2 , (5.25) 



where Wo, rn and ware the same as for qq sector (see eqs ( f4.14| ) and ( |4.15| )). After sym 



metrization over all quarks the wave-function (|5.20|) can be considered a good approximation 
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to the ground state of the A-isobar. According to eq. ( |5.25|) , the mass of the ground state 



(N — n — 1, L — I — 0) in this case is equal to M 10 io = 1.043 GeV, which is intermedi- 
ate between the masses of the nucleon and of the A-isobar. The spin-spin interaction will 
increase the mass for a total angular momentum J = 3/2 and decrease it for J = 1/2. 

In the description of Regge trajectories M ni can be considered as the effective mass of a 
diquark. According to eq. (|5.23|) for n = 1 and I = we have M\o = 0.8 GeV. This value can 



be considered as an estimation of the mass of a diquark in the ground state at n — 1 and 
/ = 0. To take into account spin-spin interaction we shall fit the baryon Regge trajectories 
considering the mass of the diquark as a free parameter. As we have already seen for mesons, 
the spin-spin interaction is particularly important in the s-wave. Therefore, if we consider 
the states with I = and L > 1, the spin-spin interaction between the third quark and the 
two others may be neglected. 

Let us classify the baryon Regge trajectories composed of u, d quarks in the diquark- 
quark picture. The lowest states are composed of s-wave diquarks: ^00(^12 = S12 = 0) and 
^11(^12 = S12 = 1). There are two nucleon Regge trajectories corresponding to the orbital 
excitations of the q — T> 00 system with / = 1/2: we denote them by N L=J _ 1 / 2 and N L=J+ i/ 2 . 
Six Regge trajectories correspond to the orbital excitation of the q — V u system. They are 
degenerated in the isospin / = 1/2, 3/2: (N - A) L=J _ 3/2 , (N - A) L=J+1/2 , (N - A)' L=J+1/2 
and (N — A)l=j + 3/2- The trajectories (N — A) and (N — A)' correspond to different total 
spins of quarks, 6*123 = 3/2 or 1/2. To describe all these trajectories, we use the two 



body equation ( |4 . 2 7| ) , where the mass of the quark was taken as for meson trajectories, i.e., 
m u = irtd = .22 GeV and the masses of "Doo and £>n diquarks have been found from fits 
to the main N and A trajectories, Nl=j_i/2 and A^ = j_3/ 2 (see Figs. 5 and 6), yielding 
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Moo — -44 GeV and M\\ = .80 GeV, which are consistent with the values given above. 

The agreement of the theoretical predictions with the available experimental data is quite 
satisfactory. We are able also to predict the A Regge trajectory composed of s — V Q0 orbital 
excitations. Our predictions (see Fig. 6) are in very good agreement with data. 

VI. CONCLUSIONS 

Starting from the BS equation in ZRA, we have introduced a null plane wave-function 
and have derived a three-dimensional two- body equation. We have generalized this equation 
to the case of finite range interactions and have discussed the choice of angular momentum 
operators in this scheme. We have considered also the three-body problem in the same 
approach, showing that in ZRA it is possible to derive a three-body null plane Hamiltonian 
with two-body interaction which satisfies cluster separability condition. We have generalized 
also this approach to the case of two- and three-body interactions with finite range. As 
an application of this scheme we have considered Regge trajectories for relativistic two- 
and three-quark systems. We have shown that linear trajectories can be derived when the 
interaction term in the mass square operator is of the oscillator form. 

We have discussed the relation between this kernel and the Hamiltonian used by Godfrey 
and Isgur J3|. Then we have applied our model to the description of the meson Regge 
trajectories q — q, q — s (s — q) and s — s. The model describes, with only four free parameters, 
m u = irid = m, m s , Wo and u, all orbitally excited meson states in a quite satisfactory 
way. We have considered also the choice of the kernel in the three-body case. It appear 
that the choice of two-body forces in the relativistic three-body equation satisfying the 
cluster separability property leads to an equation which is essentially different from the 
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semi-relativistic generalizations of the Schrodinger equation. 

We have also considered a kernel which corresponds to a three-body force with string 
junction; this operator can be taken in a form which commutes with the angular momen- 
tum operator in the free particle representation, so that the angular condition is satisfied. 
Choosing this model of interaction, we consider Regge trajectories for three-quark systems. 
We calculate all baryon Regge trajectories composed of u and d quarks in the diquark-quark 
configuration introducing phenomenologically hyperfine splitting through the difference in 
masses of the two color antitriplet diquarks, £>oo and T>\\. Our predictions are in good 
agreement with iV and A Regge trajectories. Also the prediction of the A Regge trajectory, 
which corresponds to s — V 00 orbital excitations, is in very good agreement with data. 

To summarize, the approximation proposed for the BS equation or for the four- 
dimensional three-body equation turns out to be a covariant Hamiltonian dynamics [0, 
provided we consider a convenient representation of the angular momentum operators and 
choose a kernel which commutes with such operators. Moreover, concerning quark dynam- 
ics, this approximation may be implemented by three-body interactions with string junction 
- whose importance has been already stressed in the literature - to yield a unified description 
of meson and baryon Regge trajectories. 
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FIGURE CAPTIONS 

[Fig. 1] - Meson Regge trajectories made of u- and/or <i-quarks (M[ 2 ). 

[Fig. 2] - The parameter a n \ as a function of g n i = 2(n — 1) + /, found from the solution 
of eq. ( f4.21| ) for mi = .22 GeV and m 2 ranging from 0.22 to 1.22 GeV. 



[Fig. 3] - Regge trajectories for strange mesons (M L \ 2 ). 

[Fig. 4] - Regge trajectories for ss mesons (M{J Si2 ). 

[Fig. 5] - Non-strange baryon Regge trajectories. 

[Fig. 6] - A Regge trajectory, corresponding to the orbital excitations of the s — T>„ 
system. 
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